Reading and writing charge on graphene devices 
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We use a combination of charge writing and scanning gate microscopy to map and modify the local 
charge neutrality point of graphene field-effect devices. We give a demonstration of the technique by 
writing remote charge in a thin dielectric layer over the graphene-metal interface and detecting the 
resulting shift in local charge neutrality point. We perform electrostatic simulations to characterize 
the gating effect of a realistic scanning probe tip on a graphene bilayer and find a good agreement 
with the experimental results. 



Gate tuneable conductivity underpins the operation of 
a multitude of applications envisaged for graphene, rang- 
ing from switching in mass-produced field-effect transis- 
tors to single electron manipulation in quantum dots. 
Topographic, morphological, and electrostatic variabil- 
ity often compromises device performance, so developing 
and adapting techniques for imaging and modifying lo- 
cal electronic properties continues to receive considerable 
interest. One such technique is charge writing (CW), 
which uses the biased tip of a scanning probe microscope 
to deposit charge in a dielectric layer electrostatically 
coupled to a nanodevice, mainly those fabricated from 
two-dimensional electron gas-based systems in semicon- 
ductor heterostructures [U [2] . Scanning gate microscopy 
(SGM) employs the same tip as a mobile gate for imaging 
the charge-induced changes in local conductivity. While 
SGM has also recently proven capable of making charge 
neutrality point (CNP) maps of inhomogeneously doped 
graphene flakes [3 11| , CW has yet to be used to modify 
the local CNP. In this letter we measure and modify the 
local CNP using combined CW/SGM, paving the way 
towards in situ electrostatic patterning of arbitrary po- 
tential landscapes in graphene field-effect devices. 

Our graphene devices are fabricated from flakes which 
are mechanically exfoliated from natural graphite onto a 
highly doped Si substrate capped with a 300 nm thick 
Si02 layer. The number of layers is identifled from 
their optical contrast and 5/30 nm Ti/Au contacts are 
patterned using e-beam lithography, thermal evapora- 
tion, and standard PMMA lift-off processing. For im- 
proved stability and sensitivity we operate our scanning 
probe microscope in non-contact mode (see [5] for de- 
tails) and under vacuum (10~^ mbar) with an average 
tip-surface seperation of 5-20 nm. The SGM setup used 
to image the CNP is shown in Fig. [ija). To bene- 
fit from the high signal-to-noise ratio achievable using 
a.c. detection, we modulate the voltage Vt on the tip 
(Nano World ARROW-NCPt) at low frequency (typically 
3 V @ 1 kHz) and detect the modulation of Ids (~ 250 
fiA) using a lock-in amplifier. The demodulated com- 
ponent SIds is proportional to the local transconduc- 
tance = dlos/dVT (« 0.2 iiA/V) [3l|5]. Our earlier 
work established that the back-gate voltage which nulls 
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FIG. 1: (a) Circuit used to perform nulling scanning gate mi- 
croscopy in four-terminal configuration on a graphene bilayer 
[LIA: Lock-In Amplifier, PI: Proportional- Integral controller.] 
(b) Difference between the nulling back-gate voltage and the 
bulk neutrality point voltage as a function of d.c. voltage 
applied to the tip. (c) Map of the neutrality point Vo of 
the bilayer. (d) Horizontal (dashed) and vertical (solid) line 
profiles. Arrows correspond to the n-type doping features 
indicated in the map. 



the transconductance is the local charge neutrality point 
(Vq), and spatial maps of Vq were generated by perform- 
ing point spectroscopy and extracting Vq at each point 
3J. Owing to the sign change of the transconductance 
either side of Vq, g,n{VBG) is highly suited as an error sig- 
nal which can be nulled in a feedback loop for real-time 
tracking of Vq, similar to the compensation of the con- 
tact potential in scanning Kelvin probe microscopy [5]. 
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FIG. 2: Main figure: calculated induced electron density. In- 
set: schematic of tip geometry {h — 20 nm). The tip is mod- 
elled as a blunt cone with an opening angle of 30° and a 
spherical end with R = 100 nm as shown. 

Fig. [l|a) shows the circuit used to implement "nulling" 
scanning gate microscopy (NSGM). The modulation in 
current 6Ids is fed into a software controlled feedback 
loop which adjusts the back-gate voltage to maintain a 
constant setpoint SIds{Vo)=0- Fig. [TJc) shows a map 
of Vq captured using NSGM in four-terminal configura- 
tion on a bilayer flake, revealing up to 2 V variations 
in Vq. Most notably, regions with pronounced n-type 
doping exist around the Ti/Au contacts. Line profiles 
shown in Fig. [l|d) show that the doping extends hun- 
dreds of nanometers away from the contacts. We have 
performed NSGM on other flakes in both two and four 
terminal configurations and observe similar doping pro- 
files around contacts, in good agreement with various 
microscopy studies [THHI and theoretical analyses taking 
into account the different work functions of graphene and 
titanium [TP] . 

The features in the NSGM images shown in Fig. [l|b) 
vary over hundreds of nanometers, indicating a large ef- 
fect from the conical section of the tip. This is not 
unexpected as the tip used in these experiments has a 
relatively large half-cone angle 9q w 30°, which starts 
to dominate for tip-surface separations >5 nm pT[ I12j . 
A similar long-range SGM response was reported in the 
context of SGM on carbon nanotubes [H] and subsurface 
two dimensional electron gases [131 [T3] , where it was dis- 
tinguished from finer conductance variations which occur 
on the scale of the tip apex 

In addition to the length scale of the perturbation, an 
important parameter characterizing SGM is the carrier 
density induced under the tip, which can be quantified 
through the effective capacitive coupling /3 = An/Vr, 
where An is the change in carrier density under the tip. 
NSGM enables us to measure /3 directly by sweeping the 
DC voltage applied to the tip and tracking the nulling 
back-gate voltage. Fig. [T|b) shows a plot of the differ- 
ence between the nulling voltage and the bulk neutrality 
point voltage (Vjvp) as a function of Vr. Using the known 
relationship between carrier density and back-gate volt- 
age, n = a{VBG - Vnp) (a = 7.2 x IQi" cm'^) [Tg, we 



infer a coupling between the tip and the graphene of /3w 
13 X 10^ cm^VV. 

To obtain a theoretical estimate for the length scale 
and coupling strength of the tip perturbation, we solve 
Laplace's equation numerically in cylindrical polar co- 
ordinates for a realistic tip-sample geometry. We in- 
clude the screening charge induced in the graphene self- 
consistently within the Thomas-Fermi approximation. In 
our model, a blunt tip (see inset to Fig. [2| is positioned 
20 nm above a clean, homogeneous bilayer graphene sheet 
which lies on a 300 nm thick dielectric substrate. Ex- 
perimentally we are close to the neutrality point so, in 
our idealised model, we fix the voltage at the bottom 
of the substrate to V. At the upper boundary of the 
dielectric, the electric displacement is discontinuous by 
an amount determined by the charge density induced 
in the bilayer. We calculate this self-consistently. The 
Fermi level, Ep(r), in the bilayer varies spatially with 
the voltage in the graphene and the induced charge den- 
sity, n,(r) = -fiiEpir) - Egir))/i2nP^) for Ep > Eg/2, 
is derived from the low energy parabolic dispersion for a 
bilayer [S]. Here, P = 0.539 eV nm [17], 71 = 0.39 eV 
[IS] and Eg{r) = [e^V^all i'^^'^u+lDY''^ [IH] is the band 
gap induced in the bilayer due to the bias, Vu, across it. 

The calculated charge density is shown in Fig. [2j In 
our calculation we extend the system in radial (to r = 50 
/im) and vertical (to z — 6.3 /im) directions until n.i(r) is 
converged to within 0.05% for r < 5 /im. Qualitatively, 
ni{r) is not sensitive to the band-structure model. Even 
changing the bilayer to monolayer graphene has an effect 
of less than a few percent. Instead, it is the tip shape and 
Vt that primarily determine riiir). Within r < 5 /xm we 
find that the charge density is well fitted by the sum 
of two Lorentzians. The Lorentzian peak height varies 
linearly with Vr, but the half widths are insensitive to 
tip voltage and vary by < 1% for 0.5 < Vr < 3 V. At 
Vt = 3 V the wider of the Lorentzians has a peak height, 
ni(0) = 14 X 10^ cm~^, and a half width of 3.6 /xm. This 
fits well with the observed micron-sized perturbation and 
gives a value for (3 = An/Vr = ni{Q)/VT w 5 x 10^ 
cm~^/V, similar to that found experimentally. 

We now show how the tip can be employed to mod- 
ify the CNP map by performing CW to induce holes in 
the graphene and thereby partially compensate for the n- 
type doping around the contacts. We spin coated a 100 
nm thick layer of poly(methyl methacrylate) (PMMA) to 
act as the CW dielectric over the device shown in Fig. 
[ga). Fig. [gb) shows the NSGM map in the presence 
of the PMMA layer. Note the reduced spatial resolu- 
tion due to the necessity of scanning 100 nm from 
the surface, and the enhanced n-type doping, especially 
next to the contacts where Vq reaches as low as w -9 V. 
Further work is required to use thinner and less invasive 
dielectrics. Prior to CW we acquire a high resolution 
NSGM image of Vo{x,y) next to the contacts [Fig. [3]jc), 
image 1] . To perform CW in PMMA we follow the proce- 
dure presented in |19| for creating electrostatic templates. 
The tip is brought in to hard tapping-mode contact over 
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FIG. 3: (a) Setup used to write charge in a 100 nm thick layer 
of PMMA spin-coated over the graphene. Voltage pulses (20 
Ypk-Pk @ 50 Hz) are applied to the tip in tapping mode 
to write charge, (b) NuUing scanning gate micrograph show- 
ing the charge neutrality point distribution prior to writing, 
(c) Procedure for combining NSGM with charge writing to 
sequentially read and write charge neutrality point distribu- 
tions over graphene devices, (d) Line profiles of the neutraUty 
point and frequency shift for the NSGM and EFM images, re- 
spectively. 



the doped region and a square-wave signal between -10 
and -30 V at 50 Hz is applied to the tip. These values 
bring the pulses over the « 15 V threshold pulse ampli- 
tude for charge deposition in PMMA [I^ . As expected in 
tapping mode and pulsed writing, topographic imaging 
did not show any subsequent damage to the PMMA [19] . 
For imaging the deposited charge we employ electrostatic 
force microscopy (EFM), and Kelvin probe microscopy 
when quantitative estimates of the surface potential are 
required. An EFM image showing the shift in frequency 
(A/) of the cantilever over the charged region is shown 
in Fig. [3]jc), image 2. The charge is well localized and 
corresponds to a surface potential of « -2 V. We re-map 
Vo in NSGM after CW (Fig. ^c), image 3) and construct 
the difference (Fig. ^c), image 4) to enhance the CNP 
shift. Line profiles in Fig. |3][d) show a good correspon- 
dence between the peaks of written charge and the CNP 
shift, which is w 800 mV. We found the charge to be sta- 
ble for several days under vacuum at room temperature, 
in good agreement with previous observations |19j . 

In conclusion, we have introduced a method combin- 
ing charge writing and scanning gate microscopy to read 
and write the local charge neutrality point on a graphene 
device. We mapped the variation of doping across a bi- 
layer flake and found micron-sized n-type doping from 
Ti/Au contacts. We performed electrostatic simulations 
to analyze the interaction between the biased tip and the 
graphene and deduced a capacitive coupling of /3w 5x 10^ 
cm~^/V which is in reasonable quantitative agreement 
with the value determined experimentally. We deposited 
charge with a surface potential of « 2 V and shifted 
the CNP of a region next to the metallic contacts by « 
800 mV. The use of thinner, less invasive dielectrics with 
higher permittivity would enable the in situ fabrication 
of graphene devices with arbitrary potential landscapes, 
and band-gaps in the case of graphene bilayers [20-. 
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